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Abstract 

We look at properties of vortex solutions of the extended CP^ Skyrme-Faddeev model. 
We show that only holomorphic solutions of the CP^ model are also solutions of the 
Skyrme-Faddeev model. As the total energy of these solutions is infinite these solutions 
should be interpreted as describing time dependent vortices. We describe their dynamics 
and, in partcular, point out that one of the terms in the energy density is related to the 
Noether charge of the model. 



1 Introduction 



In spite of many efforts the final description of the strong coupling sector of the Yang-Mills the- 
ory still remains a challenge for theoretical physicists. An interesting approach to this problem 
is based on some non-perturbative methods which involve soliton solutions. The importance of 
solitons in particle physics has increased significantly since it became clear that they could play 
a role as suitable normal modes in the description of strong coupling regime for some physical 
theories [IJ. The presence of solitons is an evidence of a high degree of symmetries of the 
model that can be associated with conserved quantities. The relation between symmetries and 
conserved quantities is often described in terms of the Noether theorem. However, this is not 
the case for solitons since the symmetries associated with them are neither symmetries of the 
Lagrangian nor of the related equations of motion. In fact, at the first sight there are no obvious 
symmetries that can be associated with such a large number of conserved quantities. For this 
reason they are called hidden symmetries. The study of some two dimensional integrable field 
theories has shed new light on this issue showing that hidden symmetries are, in fact, symme- 
tries of the zero curvature condition known also as the Lax-Zakharov-Shabat equation [i2|. The 
question that rises is if there exists a counterpart of such structures in higher dimensional field 
theories and, in particular, in gauge theories in (3 + 1) dimensions. An answer to such question 
was given in [3] where a generalization of some ideas of integrability to higher dimensions was 
proposed using the construction of flat connections on a loop space. Such an approach has al- 
ready led to successes in some field theoretic models and so this approach has been reformulated 
further in IH . Recently (Si [6|| this approach has been used to construct integral formulations of 
the classical equations of motion of Yang-Mills theories in (3 + 1) dimensions, Chern-Simons 
theories i n (2 + 1) dimensions, and integrable field theories in (1 + 1) dimensions. The ob- 
tained results show some interesting connections between gauge and integrable field theories. 
For these reasons this approach seems to be a useful tool for dealing with effective models for 
strongly coupled Yang-Mills theories such as the Skyrme-Faddeev model Q. 

Recently, two of us flSII, have proposed a version of the (3 + 1) dimensional Skyrme-Faddeev 
model which differs from the standard Skyrme-Faddeev model in two aspects. The first dif- 
ference is the target space of the model. In the case of the model discussed in [[8l the tar- 
get space is the complex projective space CP^ which is isomorphic to the symmetric space 
SU{N + 1)/SU{N) (g) U{1). Due to the arguments presented in m, for > 2 this could be 
an alternative description of the relevant low-energy degrees of freedom which, in the paper of 
Kondo flu, is refered to as the minimal case. The other case, for which the target space is the 
coset space SU{N + l)/f/(l)^ which is frequently used in the literature, corresponds to the 
maximal case. The version of the Skyrme -Faddeev model and its exact vortex solution in the 
simplest case of the CP^ target space has been already studied before in [10]. The second dif- 
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ference is related to the presence of some additional quartic terms in the Lagrangian. If fact, as 
shown by H. Gies in [11], such quartic terms appear unavoidably in the Wilsonian effective ac- 
tion of the SU (2) Yang-Mills theory calculated up to one loop level. This observation suggests 
a possibility of the appearance of similar quartic terms in the effective action of the physically 
more relevant case of the SU (3) Yang-Mills theory. When compared with the = 1 case, the 
Skyrme-Faddeev model on the CP^ target space {N > 1) allows us to add a further quartic 
term. For the model with > 1 this term provides a new contribution to the total energy den- 
sity in the integrable sector of the model. This extra contribution drops out when = 1. It is 
important to gain a better understanding of the main consequences of such quartic terms in the 
model, and this is one of the motivations of our study which is described in the present paper. 

In this paper we address these issues and study some exact vortex solutions of the extended 
Skyrme-Faddeev model for the target space CP^ . As shown in [[SI such theory possesses an 
integrable sector with an infinite number of local conserved currents. This integrable sector is 
defined by some specific conditions on the gradients of the fields and some relations among the 
coupling constants. We consider a large class of solutions of the sector corresponding to field 
configurations which are arbitrary functions of the variables z = x^+ix^ and ?/+ = x'^+x^, with 
x^, = 0, 1, 2, 3, being the Cartesian coordinates of the four dimensional Minkowski space- 
time. Apparently, there are no finite energy solutions inside this class, but there are physically 
interesting configurations corresponding to the vortices which are parallel to the a;^-axis with 
waves traveling along them with the speed of hght. In this paper we generalize the results of tSl 
by constructing solutions corresponding to many such vortices distributed symmetrically on the 
x^ x^-plane, and in some cases rotating relative to each other. We study some of their properties; 
such as their topological charges and the dependence of the energy per unit of length on their 
separation and on the frequency of their rotation. It is worth mentioning that our vortices are 
also solutions of the pure CP^ in (3 + 1) dimensions and so are related to the configurations 
considered in [fT2l[T4l . 

The paper is organized as follows. In the next section we briefly discuss the model, its 
integrable sector and recall some of its exact solutions and also their relation to the solutions of 
the CP^ model. The third section contains new results obtained from the study of the quartic 
terms. In the fourth section we present some examples of the exact solutions that illustrate the 
results of the previous section. We also present a numerical example for which the position of 
the minimum of the energy density depends significantly on the values of the parameters of this 
solution. The paper finished with some conclusions. 
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2 The extended Skyrme-Faddeev model 



2.1 General remarks about the model 



The model considered in dH, is an extension of the (3 + 1) -dimensional Skyrme-Faddeev model 
on the target space CP^, obtained by the inclusion of two terms which are quartic in derivatives 
of the fields. The original formulation of the model explored the fact that CP^ is a symmetric 
space, namely SU{N + 1)/SU{N) ® U{1), and so the fields where parametrized by the 'so- 
called principal' variable X (g) = ga {gY^, with g G SU{N + 1), and a being the order 
two automorphism under which the subgroup SU{N) f/(l) is invariant. Therefore, one has 
X (g h) = X (g),if h E SU{N) f/(l). For the purpose of this paper it is enough to express 
the Lagrangian in terms of the quantity P^ defined as X^^d^X = a {g) P^a {g)~^. In this 
formulation the Lagrangian of the model takes the form 



~MP') + -Jr{[P„P.]r + ^ 



TriP') +7[Tr(P^P,)]^ 



where P^j can be parametrized by the set of complex scalar fields 
lowing way: 



1 2 



P„ 



2i 



1 + • M 



(JNxN 

d^u^ ■ A 



A ■ d^u 



with 



u 



Ui 



(2.1) 

,x^) in the fol- 



(2.2) 



\Un J 

where -d = ^/l + u'f ■ u. Thus 



The symbol A denotes the hermitian matrix Ay = ?9 6ij — 
in this parametrization of the A^-dimensional complex projective space CP^ we choose to use 
complex scalar fields Ui. 



The first term in (|2.l|), namely, — ^Tr(P^P'^) corresponds to the CP^ Lagrangian which. 



in the standard formulation, takes the form 

£cPN = M\D^ZyD^Z, 
where the covariant derivative is defined as D,,^ i 



Z^ -Z 



(2.3) 



(Zt ■ d^Z)<i/. In (2.3) we have 



also included a dimensional constant since we are interested in comparing the results of the 



model based on this term alone to those based on the full model (2.1 ). This correspondence is 
established via the parametrization 



un) 



1 + luiP + . . . + \un\ 



(2.4) 



for which Z, by definition, satisfy the condition Z'^ ■ Z = \. One can check that the part of the 
Lagrangian £ proportional to is proportional \.o Lcpn i-e- ^ap = '^^cpn- 
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The last three terms in ( |2.3| ) are quartic in powers of derivatives. Amongst them the term 
proportional to 1/e^ is the standard Skyrme-Faddeev quartic term whereas those proportional 
to (3 and 7 constitute an extension of the model. The extension contains all possible quartic 
terms that can be considered for this model. Note that in the case of the CP^ target space the 
three quartic terms produce only two types of terms, namely (9^ui(9^u*)^, and {d^Uif {djju\)^ . 
Thus, in this case, out of the three coupling constants (3 and 7, only two are independent. 
This statement is not longer true for N > 2, where the interplay between the Lorentz indices 
IX and u, and the internal indices i and j produces extra independent terms. Hence, in the case 
N >2 both /3 and 7 can play a significant role. 

One of the most important results presented in ^ was the demonstration that the extended 



Skyrme-Faddeev model (2.1) possesses an integrable sub-sector (sub-model) with an infinite 
number of conserved currents. The integrability condition has the form 



d^Uid^Uj = 0, for any i, j = 1,2, . . . , N. 



(2.5) 



In what follows we shall refer to (|2.5[) as the constraints. As shown in ^ the equations of 



motion of (2.1 ) together with the constraints ( 2.5[ ) imply that for any arbitrary functional G of 



the fields Ui and u*, but not of their derivatives, there exists a current which is conserved. 



i.e. d^'J^ = (see M for details). 



In addition, it was also shown in [[8l that if one imposes a particular relation among the 
coupling constants; namely: /3e^ + 76^ = 2, the integrable sub-model possesses a wide classes 
of exact solutions. These solutions are given by arbitrary meromorphic functions of the complex 
variable z = + ix^ and of the real variable ?/+ = + x", with x^, /i = 0, 1, 2, 3, being the 
Cartesian coordinates of the Minkowski space-time. Indeed, the configurations 



Ui = Ui[z,y+} 



Ui[z,y+} 



(3e^ + 76^ 



(2.6) 



satisfy the Euler -Lagrange equations corresponding to (2.1 ) as well as the constraints (2.5). In 



fact, there is more to it. The configurations (2.6) satisfy the following equations 



d'^Uid^Uj = 0, d^[d^Uidf,Uj] = 0. 



(2.7) 



Note that solutions of (2.6) are also solutions of the CP model in (3 + 1) dimensions, i.e. 



the model defined by the Lagrangian (2.3 ) or the first term in (2.1 ). That fact has been explored 
in [fT2l| to construct vortex solutions for the CP^ model. Therefore, a key point to be stressed is 
the observation that the CP'^ model and the integrable sector of the extended Skyrme-Faddeev 



model share a wide class of solutions given by the solutions of (2.7 ). The vortex-type solutions 
presented in [|T2l are more general that those of [8]. Moreover, the vortices can be located at 
any distance from each other - and so we can study their dependence on this distance (in the 
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previous case this distance was zero as they were located on 'top of each other'). In fact we find 
that the energy per unit length of the vortex depends in a nontrivial way on the distance between 
any two or more individual vortices. We would like to add, however, that not all solutions of 
the CP^ model are simultaneously solutions of the extended Skyrme-Faddeev model. The 
class of the so-called mixed solutions, Ui{z,z,y^) of the CP^ model, see also [fT3ll . that we 
have discussed in [|T4ll . does not satisfy the constraint (|2.5[) and so they are not solutions of the 



integrable sub-model. In this paper we concentrate our attention on the contribution to the total 
energy density generated by the quartic terms. We are especially interested in the dependence 
of the energy per unit length on the distance between vortices. 

2.2 The Hamiltonian 

It is convenient to split the Hamiltonian into three parts 

= sM^n^'^ + -H^^)) + 64(7 - /3)-h(=^), 

where the subscript c implies that the Hamiltonian is already restricted to the integrable sector 
by assuming ( 2.6| ). The first two contributions are given by 



(1 + Mt-M)2 ' (l + Mt.M)2 ' ^^-"^ 

where Af^ = (1 + • u)5ij — UiU*. The quartic contribution can be cast in the form 

1 ^ 

where Bji = {dzUjdy^ui — d^uidy^Uj) and 5*^. = {dzU*dy^u% — dzuldy^u*). The reduced 
Hamiltonian "He is positive definite for 7 > /3 since each of its contributions "H*^"), a = 1, 2, 3 is 
positive definite. The proof of the positive definiteness is given in ||8l . 



The Lagrangian (2. 1 ) is invariant under the phase transformations Ui — > e^^^Ui, i = 1, . . . , N. 



The associated Noether currents (for the restricted sub-model) are given by 

where we have split the expression of J^*^*) into two parts 

Jll"^ - Y^ il + uUuy f^ [<(A^)..5.^. - d,u*{A%u;\ (2.10) 
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and 



1 1 ^ 



- {A'),,{A%kUid''uk{d,u*d,u*i-d,u*d,u*)]. (2.11) 

Next we observe that both contributions J^*) and jTjf*-' are conserved independently as a 
consequence of the fact that in the integrable sector the fields satisfy the relations (2.7 ). Indeed, 
using ( |2.7 ) one can check that d^J^^^ = 0. Since the Noether theorem implies d^J'^^^^ = 
it follows that d'^J'^''^ = 0. Another way of obtaining this result is to observe that (2.7) also 
defines an integrable submodel of the pure CP^ model, and that jrj*^ is the Noether current 
of such a model, associated with the symmetry Ui — )■ e'"^'Ui, i = 1,. . . ,N. Therefore, J'lf'^ is 
conserved, and so is Jj^) by the same argument. 

Following [8J we shall express the restricted Hamiltonian density Tic in terms of the time 
component of above mentioned Noether currents. We find that 

^0 = ^- n\ t %2 (2.12) 



and 



1 1 ^ 



Ji'^ = W^jT—i — u E [i^^)iji^'^)kiu*d,ul{d+Ujd,ui-d,Ujd+ui) 

-\-U ■ U) j^^^i^i 

- {A^),,{A%kU^d,Uk{d+uW,u*,-d,uW+ul)]. (2.13) 



Thus the Hamiltonian density can be rewritten as 

N 

= SM^dA In (l + u^.u) + 8M2 J2 hJ^'^ + 64(7 - /3) ^ hJc 



N N 

5 

i=l i=l 



where the first term is purely topological and the last two, related to the Noether currents, 
involve derivatives w.r.t. y^. Thus it follows that the energy per unit length can be cast in the 
form 

AT 

£ = SttM^Qtop + E (SttM^qW + 647r(7 - /3)Q») , 

i=l 

where the expression in the bracket describes the Noether charges Q^''\ The topological charge 
is given by the integral 

Qtop = - f dx^dx^n^^^ (2.14) 



and in a similar manner we have introduced Q^^^ and Q*^*) i.e. 

Q« = - [ dx'dx^H^'\ Q« = - / dxVU^^^ (2.15) 
The reason why we split the charges Q^^^ into two parts will be made clear in the next section. 
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3 The relation between the Noether charges and the topolo- 
gical charge 



In this section we restrict our study to the case N = 2. The reason for this is twofold. Firstly, 
this case is more interesting from the physical point of view since the Yang-Mills theory is a 
gauge theory based on the symmetry group SU{3). Secondly, it is the simplest case for which 
Ti^^^ is not identically zero. 

In contradistinction to the parts H^^^ and 'H^'^\ the contribution to the energy density H^^^ 
contains terms that have both type of derivatives i.e. w.r.t. z and The term is a total 
derivative and therefore is purely topological whereas V}-'^^ has nothing to do with topology 
since it does not even contain any derivatives w.r.t. z or z. From the point of view of topology 
one can, however, expect that 'H'-'^-' may have some interesting properties. Indeed, encouraged 
by the results of our numerical studies we have managed to establish a relation between the 
topological charge Qtop and the contribution Q^*^ to the total Noether charge. The topological 
charge arises in the contribution to the energy per unit length calculated from "H*^^) when only 
one of the functions ui or u-2 depends on z. Without any loss of generality we can choose, for 
instance, Ui = ui{z, ?/+) and U2 = U2{y+). One can then check that in this case the expression 



, given in \2.9\ , takes the form 



d2u\d,ui 
[1 + w ■ uY 



(3.16) 



P, we observe that f.^^^f"'". 



(o+|mi|2)3 



1 ^ 

2 da 



1 dzvldzVi 
a {l+|i;i|2)2 



Denoting a = 1+ | U2 
where we have denoted vi = ui/y/a. When integrating over and we can use the fact 
that the integral dx^dx'^-^^^^l^:^^^ is invariant under rescaling of Vi by a constant, since it is a 
topological quantity. Therefore, we obtain the result 



/ dx^dx^n^^^ = vrgTop- 

Ji?2 I 



d+uld+U2 

1+ |M2p)2' 



(3.17) 



where Qtop is defined in (2.14). 

Thus we see that the energy per unit length for this special choice of scalar fields ui and U2 
is a product of the topological charge and of a function which depends only on y^. This implies 
that, for some special solutions, the topology can play an important role also for the energy of 
the interaction generated by the quartic term. 



The zero-components of the Noether currents (2.13 1 in this case become 



(1) 



0, X 



(2) 



{1 + ■ uy i 



[uld+U2 - 



(3.18) 
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The and integration then gives 

/ d^^da^jf ^ .Q(" = .Qr.,^ '"'^';": : (3.19) 

Jij2 2l (1 + \U2\ ) 

which is, again, a product of the topological charge and of some function of the variable y+. One 
can easily check that the only dependence on of U2{y+) that guarantees the proportionality 
of the two last integrals is U2{y+) = ae''^^+ where a and A are real constants. In order to see 
this one can express the complex scalar field U2 in terms of two real fields R{y+) and ^{y+) i.e. 
U2 = R{y+)e''^^y+\ Then it follows from the proportionality condition that 

d+U2d+U2 1 uld+U2 — U2d+U*2 

A- 



and so that 

R'^ + RV = Ai?2$'. 
This equation is satisfied by R being a constant R = a and $ = Xy+. Then 

Q^"^ = Qto,jj^^,. (3.20) 
(1 + a"^)^ 

It is easy to observe that this is the only case when the ratio of the two charges is a constant. 
In consequence, the energy contribution can be expressed in terms of the corresponding Noether 
charge. If one choses, for instance, U2{y+) as a real function U2 = R{y+) then the contribution 
to the energy per unit length remains some function of y^, while still keeping its proportionality 
to Qtop, whereas both Noether charges vanish. 

It is clear that the other possible choice, i.e. ui{y+) and U2{z, y+), is equivalent to the one 
discussed above. The third possibility leading to the proportionality of the energy per unit length 
/ dx^dx'^'H'^'^^ to the topological charge Qxop arises when both scalar fields ui and U2 possess 
the same dependence on z (or z) i.e. Ui{z, y^) = v{z)gi{y^). This can be easily seen if, instead 
of parametrization in terms of ui and U2, one considers Z, introduced in ( 2.3[ ), parametrized as 



^ (/^l; /^2) /"s) (3 21) 



Comparing this with (2.4) we note that there are several equally good possibilities of the choice 



of Ui e.g. Ui = fi/ fs or ui = /1//2, U2 = fi,/ f2 e.t.c. This demonstrates that two possibilities 

. , _ N / U2 1 

(mi,U2,1) and (1, ?i2, mi) = (1, — , — ) 

Ui Ui 



are totally equivalent. Moreover, this discussion shows that in this new reparametrisation we 

, ,\ , andM2(^) = 

'"\^)9i{y+) ^ ' 9i{y+) 



are taken back to the previous case Ui{z,y+) = , , — ^ and U2iz) = Hence these two 



results are equivalent. 
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4 Some examples 



In this section we give some examples of solutions that are based on the solution of the CP"^ 
model discussed in [fT4ll and some generalizations of the multi vortex solutions discussed in 
BHl. We concentrate on the topological case (the quartic term becomes a topological one) but in 
the final subsection we discuss also a numerical solution which is not of this type. For such a 
solution the minimum of the energy per unit length depends on the coupling constants. 



4.1 Example 1 



For some arbitrary choices of the set of functions fi{z, y^) introduced in (3.21 ), or equivalently 
Ui{z, y+), the energy density per unit length is a very complicated function. Therefore it is quite 
hard to discover what is the relation between the values of the parameters of the solution and 
the distances between maxima of the energy density. In order to overcome this problem it is 
convenient to study some cases with a symmetric distribution of maxima of the energy density. 
The positions of maxima depend on zeros of the functions fi{z, ?/+). When only one of these 
functions has zeroes then maxima are located at these zeroes. For this reason we first consider 
the following case 

f2{z,y+) = a^z + a^e^^y^ 

h{z,y+) = a4, (4.22) 

where N <Q. 

The power Q is the maximal degree of polynomials in z which so gives Qtop = Q- When 
02 = the energy per unit length / dx^dx'^'H^^^ has to be proportional to the topological charge. 
Moreover, in this case U2 = ^^6**^^^+. According to the previous section the Noether charge and 
the topological charge are proportional for 02 = 0. The energy density contributions thus take 
the form 

H« = |^, n^'^ = ^^, H(^) = |^, (4.23) 

^iQ,N) ^(Q,N) ^(Q,N) 



where 



X(Q,^) = al + al + ay^ + r^'^ + 2a,r''+'^ COS [hy+-{Q~N)^], 



+ 2NQair^^^-^ cos [kiy+ - {Q - N)^] 
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{Q,N) 



W(^Q^N) = alalkj + a\[a\kl + a^(fci - k2f]r'^^ + 03^2^^'^ 
- 2aiask2{ki - A;2)r^+« cos [kiy+ - {Q - N)ip]. 
There are two special cases when the formulas simplify a lot; namely: 

1. N = 0, 

2. N = Q. 



4.1.1 The case N = 

For = and ai 7^ 0, each contribution T-L^^^ 'H^'^^ and T-L^^^ have exactly Q symmetrically dis- 
tributed maxima lying on circles with radii r^^\ r*^^) and r^^^ respectively. The angular position 
of k-th maximum is given by 

1 



My+) = Q [kiy+ + [2{k - 1) + ^(ai)]7r] 



(4.24) 



where 9{ai) is the Heaviside step function and k = 1,2, . . . ,Q. The radii r^'^\ a = 1,2, 3, can 
be calculated exactly. We present here only r^^^ and r^'^^ since r*^^^ is given by a very complicated 
expression. They take the form 



.(1) 



,(3) 



l«l| + VQV + (Q'-l)(«3 + «4) 



Q + 1 



(g + 2)|ai| + ^9Q^al + 4(Q - 1)(2Q + l)(ai + a|) 
2(2Q + 1) 



From these formulas it is clear that Oi is a crucial parameter which determines the distance 



between individual maxima. The leading behaviour for |ai| ^ 1 is r 



(1) 



^/'^ + ... and 



,(3) 



ai 



p-/'^ + . . . whereas r^^^ 



c\ai 



l/Q 



+ . . . where c depends on all other parameters 



ki, k2, as, 04 and Q. An example of the energy density contributions is sketched in Fig{T] For 
ai = the energy density contributions form a symmetric crater or a peak. When oi increases 
the picture distorts itself into three gradually emerging maxima which for ai ^ 1 became three 
well localized peaks. These peaks rotate with the speed of light. When visualized in 3-dim. 
space the positions of maxima of the energy density take the form of three rotating spirals. The 
calculation of the energy per unit length gives 



/ dx^dx^n^^^ = Q, -I dx^dx^n^^^ = Q 



71 JR'^ 



k2Q 



(2) 



(4.25) 
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where Qtop = Q and Q^^^ is the Noether charge. The last formula can be obtained directly 
from (3.20) for X = k2 and a = a^/a^. 

The quartic contribution vanishes in three cases: 03 = 0, 04 = and k2 = 0. For each 
of those cases there is effectively only one Ui that matters, which means that the target space 
reduces from CP^ to CP^. It is important to stress that both integrals (4.25) do not depend 
on ai and therefore also on the distance between vortices. The only dependence comes from 
the integral / dx^dx'^'H^'^\ Unfortunately we cannot calculate this integral analytically. The 
numerical integration, see Figj2[ shows that the energy per unit length takes a minimal value for 
some non-zero a^**". When 03—7-0 the position of the minimum also tends to zero, a™" — )■ 0. 

The approximate dependence on ai can be calculated analytically around ai = by inte- 
grating coefficients of its Taylor expansion 



/ dx'dx^n^^^ = ^X„af. 

n=0 



Note that only even coefficients X„ are non-zero. The results of integration give the following 
general expression 

, 2 



-1 



\n+l 



n 



X 



n+l 
.2\2 



ri 



ri 



X 



kin^Q\al + alY + klal{nQ - l)[{nQ - l)a^ - Qa 



2kik2a\nQ{nQ — l){a\ + a\ 



(4.26) 



where 



(a)„+i = a(l + a)(2 + a) . . .{n + a). 

is the Pochhammer symbol. The expansion converges for |ai| < 1. Unfortunately the analytical 
curve does not reach the minimum; nevertheless in the region of convergence it serves as a test 
confirming our numerical computations. 



4.1.2 The case N = Q 



In the case N = Q the energy density of all contributions 'H^"'\ a = 1,2,3, has an axial 
symmetry. In this case one has to exclude the values ai = ±1 since they lead to the vanishing 
of /i for values of y+ such that 6*^^^^+ = ±1. The integrals corresponding to a = 1 and a = 3 
remain unchanged and they are given by formulas ( 4.25| ). Note that the corresponding energy 



densities (4.23 1 depend on y+ via the periodic function cos {kiy^ 



It is important to stress that the formulas ( 4.25[ ) are the same in both cases = and 
N = Q but the energy density functions and their dynamics are completely different. In the 
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first case there are rotating peaks whereas in the second one there are oscillating symmetric 
rings. The only dependence of the energy density on ai comes from the integral / dx^ dx^'H^'^^ 
which also depends on y^. The integral can be computed analytically giving 



dx^dx^n^^^ = T, 



IT JR'^ 



where 



V n=l 



and 



1 - 




(4.27) 



ag = al{al + Qa1)kl + alalkf + afal^ki — 
+ a\alal{2kl - 2kik2 + Qkl) + 
— 2aia\k2 a\{ki — (5^2) + '^3(^1 — ^2) cos {kiy^ 



= (03 + 04) 



al + 2ai cos (kiy^ 



The integral (4.27) is in excellent concordance with its numerical estimate. Some corre- 
sponding curves are presented in Figj3j This integral is infinite for oi = 1 and kiy^ = (2?2 + l)7r 
or ai = —1 and kiy+ = 2nn where n = 0, ±1, ±2, . . . Since we have excluded ai = ±1 from 
the space of parameters the integral is always finite but it could take very high values close to 
the points where fi vanishes. The energy per unit lengths always tends to zero for oi — )■ ±00 
where the leading behaviour is ~ |ai|^^/*^. This example shows that for some energy density 
contributions, like a quartic term which we discussed here, the integral over the x^x'^ plane can 
be time independent in spite of the fact that the energy density is not. It is important to stress 
that there is no rotation here which explains the time independence of the energy per unit length. 



4.2 Example 2 

Another interesting solution can be obtained by a straightforward generalization of the vortex 
solution presented in [[8]|. Let us consider, for instance, the one parameter solution of the form 

Uj{z,y+) = {z- 6p{z + 5)^^e'^^y+, J = 1,2 (4.28) 

where 5 is some real number. For the values of 5 being sufficiently large the energy density 
is localized around the points 2; = ±5 taking the form of either peaks or craters. As a vortex 
is not a point-like object in order to study the dependence of the energy per unit length on 
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the separation between vortices one has to define what the expression "the distance between 
vortices" means. In order to do so one needs to chose some characteristic parameters of the 
solution. In the case of the solution (4.28) we only have the points z = ±5 in the complex 
plane. For 5^1 the value 25 is a distance between centers of two local peaks (or craters) in the 
plot of the energy density. However, for smalles 5's the picture of the energy density becomes 
more complicated and the reference to the points z = ±5 as to the centers of the vortices is 
no longer valid. Thus, we can think of 26 as the distance between vortices only for separated 
vortices. 



4.2.1 Topological sector 

According to discussion in the previous section, for the solutions of the form Ui{z,y^) = 
f (z)e*'^'^+ the quartic term contribution to the energy density has a topological nature. Tak- 
ing ni = n2 = n and mi = m2 = m one can replace the solution by the equivalent one 

which gives 

JR-^ 4 
where for 5 > 0, Qtop = + m| for nm > and Qtop = max(|n|, \m\) for nm < 0. 
When m = then Qtop = I'^l and vice versa. The integral does not depend on 5 which is 
a counterpart of the parameter ai from Example 1. Fig|4] shows the integral ^ Jj^2 dx^dx'^l-L^'^^ 
which is the only contribution to the total energy per unit length which depends on the distance 
between vortices. For solutions of the form ( |4.28[ ) not all combinations Ui and mj lead to the 



finite energy per unit length. In fact a grid of points corresponding to acceptable combinations 
of exponents forms a quite complex structure but we shall not discuss this problem here because 
it is beyond the scope of present paper. 



4.2.2 Outside the topological sector 

The common property of all solution having the quartic term topological is that the minimum 
of the total energy density is determined only by properties of the integral ^ Jj^2 dx^ dx'^'hP'\ 
Because of this the position of the minimum does not depend on the coupling constants. When 
the quartic term cannot be reduced to the topological case the minimum (if it exists) depends on 
the mutual relation between quadratic and quartic terms. The position of such a minimum is a 
function of the ratio 8(7 — /3) /M^. This possibility is new since the quartic term for the model 
with a target space being CP^ is identically zero if one imposes the integrability conditions 
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(2.5 1. As an example we consider a solution of the modified Skyrme- Fadeev model in a CP^ 



target space. Such a solution corresponds to (4.28) for the choice rii = n2 = —1, rrii = — 1, 



1712 = S, ki = 1 and ^2 = 2. FigJS] shows two contributions to the energy density. The quadratic 
contribution leeds to energy per unit length that increases with a distance between two multi- 
vortices. This contribution has a minimum at 5 = 0. The contribution to the energy per unit 
length given by the quartic term has maximum for 5 = and it decreases as 5 tends to infinity 
what means that the quartic term leads to a repulsive interaction between two multi- vortices. 
Fig shows the sum of quadratic "H*^^^ and V.'^^^ contributions to the total energy per unit length 



ttJr^ V M2 



In both cases the combination of the terms leads to a minimum situated at some 5 > 0. The 
left picture corresponds to the choice (7 — (5)/M'^ = 1/8 whereas the right picture shows the 
combination with (7 — (5) /M^ = 1. 



5 Concluding remarks 

In this paper we have shown that holomorphic solutions of the CP^ model are also solutions 
of the extended CP^ Skyrme-Faddeev model. In fact, these are very special solutions since 
they belong to an integrable sub-model of the CP^ Skyrme-Faddeev theory, defined by the 



constraints (2.5), and the condition on the coupling constants given by /3e + 76 = 2. The 



constraints ( 2.5| ) imply that such sub-model possesses an infinite number of conserved currents. 



As the total energy of these solutions is infinite these solutions should be interpreted as describ- 
ing time dependent vortices, and they generalize the results obtained in ll8l [T2l[T4l . 

In this paper we have studied the properties, and in particular the dynamics of these multi- 
vortex solutions of the CP^ Skyrme-Faddeev model. The structure of the solutions is very 
complex and diverse, and we have in fact considered only some types of solutions which are 
interesting from the physical point of view. We have shown for instance that in some cases one 
of the terms in the energy density is related to the Noether charge of the model, and in some 
other cases the energy density can be factorized in the product of two terms where one of them 
is the topological charge. 

It has been put forward in [|9| that the CP^ Skyrme-Faddeev model might describe some 
aspects of the low energy (strong coupling) regime of the pure SU (2) Yang-Mills theory. If that 
is indeed the case, the solutions constructed in this paper certainly must play a role, describing 
some type of low energy excitations of the Yang-Mills theory. For those reasons more research 
work is needed to get a better understanding of these interesting phenomena. 
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Figure 1: The Q = 3 solution with 02 = 0. The functions T-L^^^ (first column), T-L^"^^ (second 
column), T-L^^^ (third column) for oi = (first row), ai = 0.2 (second row), ai = 1 (third row), 
ai = 5.0 (fourth row) and 03 = 1.0, 04 = 1.0, ki = 1.0, k2 = 2.0, y+ = 0. 
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Figure 2: The case 02 = 0. The integral / dx^ dx^VP^ as the function of ai for 03 = 1, 04 = 0.8 
ki = 1.0, k2 = 2.0. The numerical integration are given by points whereas the curve shows a 
few first terms of expansion at ai = 0. The expansion converges for |ai| < 1. 
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Figure 3: The integral ^ Jj^2 dx^dx'^'H^'^^ as the function of time x^. The left picture contains 
some numerical values (the points) and the analytical curve. It corresponds to the choice Q = 2, 
ai = 3, 03 = 1, 04 = 3.5, ki = 1 and /c2 = 2. The right picture was obtained for the same 
values of the constants and various Q; from the top: Q = 2, Q = 3, Q = A, Q = 5. 
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Figure 4: The integral ^ dx^ dx'^'H^'^^ as the function of 5 for ni = n2 = —3, nii = m2 = 1, 
ki = 1 and k2 = 2. This is the only integral which depends on 6 since the quartic contribution 
is topological for this case. 





Figure 5: The integral ^ dx^dx'^VP'^ (left) and \ jp,2 dx^dx'^V.^^^ (right) as the function of 5 
for Til = n2 = —I, Till = —1, m2 = 3, ki = 1 and /c2 = 2. The quadratic contribution to the 
energy density has minimum at 5 = 0. 
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Figure 6: The integral ^ dx^dx^n^^^ + "H^^)) (left) and ^ dx^dx\'H^'^^ + 87/(3)) (right) 
as the function of 5 for ni = n2 = —1, mi = —1, m2 = 3, fci = 1 and /c2 = 2. The minimum 
^min > appears as the effect of mutual interaction of the quadratic and the quartic term. 
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